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We derive rigorous lower bounds for the combinations of light quark masses rris + niu and + niu- 



1. Introduction 

The masses of the light u, d and s quarks are 
fundamental parameters of the Standard Model 
and play an important role in low energy hadron 
phenomenology as well as in our understanding 
of both chiral symmetry breaking and the ori- 
gin of mass. Until perhaps a year ago, a cer- 
tain concensus had been reached as QCD sum- 
rules 1^ and quenched lattice QCD Q yielded 
compatible values. However, recent analyses of 
the cutoff dependence of quenched lattice results 
have yielded significantly lower values, with an 
indication that unquenched values could be even 
smaller [Q,^. Perhaps even more surprising are 
the results of the dynamical Wilson fermion cal- 
culation by SESAM § who find very small u and 
d masses but a strange quark mass in agreement 
with the QCD sumrule determinations. (For even 
more recent developments, see Q.) In light of this 
rather confusing situation, we have sought to de- 
rive rigorous lower bounds for the combinations 
of light quark masses iris+mu and md + mu below 
which lattice results cannot fall without implying 
a surprising breakdown of perturbative QCD. 

In our analysis, we consider two channels. The 
first, corresponding to the two-point function 

V'5(g') = * j d^xe'9-(0|r{A(x)At(0)}|0) , (i) 

with A = dfj,A'^~i{ms + mu)'.s"f'^u: or i{md + 
niu)-dj^u:, is a pseudoscalar channel. The sec- 
ond is a scalar-isoscalar channel and is associ- 
ated with 

= 1 [ d'x e'^-(0|r {5(a:)^t(o)} |o) , (2) 
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where S—ifiy.uu: ~\- idd:] and m~{mu + md)/2. 

Both two-point functions satisfy a twice sub- 
tracted dispersion relation, 

Im-0(5)(t) 



dt 



^2^3 



(3) 



/o {t - r 

which relates an integral of the hadronic spectral 
functions, Im7/'(5), to the QCD behavior of the 
two-point functions in the deep euclidean region. 
For — ^2_n2^ A 2 



Q >Aqcd. the LHS has an OPE 
1 ,„ 1 
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combination, the 
I, and V('5)2, V'(5)4 



where m is the relevant mass 
^(5)0 ^""^ known to four loops 
contain mass and condensate corrections 

The RHS is determined by the intermediate 
hadronic states which can contribute to the spec- 
tral ftmctions through unitarity: 



llmvl/(5)(g2) = (2^)3 



x(o|j(o)|r)(r|jt(o)|o) , 



(5) 



with J=A, S. For the pseudoscalar case, the 
lowest-mass contribution is |r) = |i4r+} (or |7r+) 
for the dj^u current) while for the scalar case it 
is Itttt). The K'^^n'^) contribution is determined 
by fK{Tr) and Mx(7r). The Itttt) contribution is de- 
termined by a single form factor F{t), which has 
been well studied in the low-energy region . 

It is clear from Eq. (^ that the spectral func- 
tions are sums of positive terms so that keeping 
only the lowest-mass contributions gives lower 
bounds to these functions. Then, Eq. (||) and the 
fact that the ^'"5) are proportional to the square 
of the quark masses of interest, can be used to 
translate this bound into a lower bound on these 



masses. 
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2. Pseudoscalar Channel 

Following [jll] we consider moments of the 
hadronic continuum integral (fo=(MK,^+2A/,r)^): 

— X+(7r+) contribution) . (6) 

Positivity of the spectral function then implies, 
for the first three moments, 

Sat > 0, Sm > Sm+1, > (Si)' , (7) 

with iV=0, 1,2 and Af^O, 1. 

The first bound we consider is the one aris- 
ing from Eo>0. We find, for the running quark 
masses: 

K(g2) + m„(Q2)]'> M-!H&^x (8) 

where is the scale at which the OPE of Eq. (^) 
is performed. Using three-loop running, we shall 
give all numerical results in the MS scheme at 
a reference scale of 4 GeV^. Clearly, this bound 
falls off quickly with Q^, so one would like to work 
with as small a as possible. However at small 
Q^, the OPE of Eq. (^), and in particular the 
pQCD expansion, breaks down. We consider that 
values of Q > 1.4 GeV are safe. Our results are 
shown in Fig. |l|. For low Q, the preliminary, cen- 
tral, unquenched results of |^J^ are excluded. It 
should be emphasized that in Fig. |^ we are com- 
paring a "calculation" of the quark masses with 
a bound which is saturated in the extreme case 
where all contributions to the spectral function 
other than are neglected. 

The bounds arising from Si 2>0 give an im- 
provement due to the fact that higher and higher 
moments increase the weight of the low-energy 
part of the hadronic spectral function. Some of 
this improvement is offset by the fact that the 
OPE converges less and less well so that one 
has to go out to larger values of Q. A simi- 
lar phenomenon occurs when one considers the 
quadratic bound given by SoS2>(Si)^. Never- 
theless, the net improvement is still significant, 
as seen in Fig. |l|, where we believe that values 
of Q > 2 GeV are safe. We find that even the 




1.5 2 2.5 3 3.5 4 

Q(GeV) 

Figure 1. Bounds from Eo>0 (dashed) and quadratic 
bounds (solid) for AliL=290,380 MeV. The dot- 
dashed curve illustrates the uncertainty due to un- 
known 0{ai) terms (A1£L=380 MeV). Values of 
rus +mu below the various bounds are excluded. The 
horizontal lines correspond to the central quenched 
(upper) and unquenched (lower) lattice values of 

central quenched values of |^,|^ may be excluded. 
This bound allows less trivial continuum spectral 
functions than the bound of Eq. ^: a (5-function 
with arbitrary position and weight saturates it. 

A quadratic bound can also be derived for m 
and is shown in Fig. |^. Again, the lattice results 
of l^,^!^ are in serious difficulty. Note that this 
bound is significantly more stringent than those 
of where different conclusions were drawn. 

The known radiative corrections to the Sjv are 
large (more so as N increases) and positive, and 
exhibit near geometric growth. To give some 
measure of the uncertainty associated with po- 
tentially large, higher-order terms, we also show 
in Figs. |] and H the quadratic bounds obtained by 
assuming that the O^a*) terms are equal to the 
0{al) terms to the power 4/3. Should this be- 
havior persist to even higher orders, the bounds 
would be lowered further. 

3. Scalar-Isoscalar Channel 

As noted above, the lowest mass contribution 
to the scalar-isoscalar channel is determined by 
a form factor, F{t). xPT to two loops and a 
dispersive analysis yield |ic|] : 

F{t) = F{Q)[l + {{r^):/&)t + 0{e)] (9) 
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Quadratic bounds (solid) for A^iL=290 
The dot-dashed curve iUustrates the 
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Figure 2. 
380 MeV. 

certainty due to unknown C'(Qff) terms (7 
MeV). The long-dashed and short-dashed curves cor- 
respond, respectively, to the scalar bounds with and 
without information on F's phase. The 2 curves in 
each set reflect variations of all input parameters, in- 
cluding a1£L, in their error ranges. Values of rh below 
the various bounds are excluded. The horizontal lines 
correspond to the central quenched (upper) and un- 
quenched (lower) lattice values of Q. 

with F(0)=A/2(0.99±0.02)-t-O(M6) and (r^)^ = 
(0.59±0.02) fm^ Following the methods of 
we construct bounds on m using this information. 
The QCD behavior of ip" is similar to that of 
i/'s up to very small ruu^d mass corrections which 
we neglect. The corresponding OPE appears well 
behaved for Q > 1.4 GeV. The bound we obtain 



> 



M° 12x(0.99±0.02)^ 



1+- 



(10) 



where .7-" is a nontrivial function of its arguments. 
Though this bound is quite similar to the axial 
current bound of Eq. (||), it is 1/A^c-suppressed 
and chirally suppressed by an additional factor of 
M^. One might expect, therefore, that it would 
be much worse; yet it is surprisingly competitive 
(see Fig. ||). This is due to the fact that the ab- 
sence of a prominant scalar-isoscalar resonance is 
compensated by rather large chiral loops and the 



fact that, numerically, Af^ and are not so dif- 
ferent. 

This bound can be improved by taking into 
account the phase of F{t) which is given by 
the I~J—0, TT-TT phase shift in the elastic re- 
gion AM^>t>lQM^. Incorporation of this in- 
formation requires solving two integral equa- 
tions iQ. Restricting to a conservative re- 
gion, 4Af^>t>(500 MeV)2, we use the phase-shift 
given by a resummation of chiral logs calculated 
at the two-loop level, which agrees well with 
available experimental information | |To| . The im- 
proved bound we obtain is shown in Fig. ^j. Again, 
we find that the lattice determinations of 
are in difficulty with this bound. We should add 
that the sumrule results of are all compatible 
with our lower bounds. Also, since this talk was 
given, two other attempts at putting bounds on 
the light-quark masses have been reported |jl^,|l6j . 
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